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Abstract 

The order a 2 contributions to the coefficient functions corresponding to the asym- 
metric fragmentation function F^(x, Q 2 ) in e + e~ annihilation are calculated. From 
this calculation we infer that the order (a s /4n) 2 correction to the flavour asymmetry 
sum rule is non vanishing and amounts to — 12/5oC*fC(3)- We also study the effect of 
the higher order QCD corrections on F^(x, Q 2 ) and compare them with the OPAL 
data. The latter put a strong constraint on the valence part of the fragmentation 
densities D^(x,fi 2 ). 



The measurement of the fragmentation functions in the process 

e+e- ->7,^-» H+ "X", (1) 

provides us in addition to other experiments like deep inelastic lepton-hadron scat- 
tering, with a new test of scaling violation as predicted by perturbative quantum 
chromodynamics (QCD). Here "X" denotes any inclusive hadronic final state and H 
represents either a specific charged outgoing hadron or a sum over all charged hadron 
species. This process has been studied over a wide range of energies of many different 
e + e~-colliders. Data haven been collected from TASSO [|TJ (y/s = 22,35,45 GeV), 
MARK II @ and TPC/2 7 = 29 GeV), CELLO § (^i = 35 GeV), AMY § 
(y/s = 55 GeV) and DELPHI §, ALEPH §, OPAL § (^/i = 91.2 GeV). 
Following the notation in |J the unpolarized differential cross section of process (|l]) 
is given by 

d 2 a H 3,, 2nydaS 3 . 2 n da¥ 3 do 1 ! 

1 + cos 2 6>)-5 + - sin 2 0-£ + - cos 0-A (2) 



* = =f, 0<x<l, (3) 



dxd cos 9 8 rfx 4 <ix I <ix 

where x stands for the Bj0rken scaling variable 
2pq 

Q 

and p and q (q 2 = Q 2 > 0) are the four-momenta of the produced particle H and the 
virtual vector boson (7, Z) respectively. The variable 6 denotes the angle of emission 
of particle H with respect to the electron beam direction in the center of mass (CM) 
frame. The transverse, longitudinal and asymmetric cross section in (Q) are defined 
by cr T , o L and a a respectively. The latter only shows up if the intermediate vector 
boson is given by the Z-boson and is absent in purely electromagnetic annihilation. 
Before the advent of LEP1 the CM energies were so low (y/s < M%) that a a could 
not be measured and no effort was made to separate from a? so that only data 
for d 2 a H I dxd cos 9 (0) were available. Recently, after LEP1 came into operation. 
ALEPH [0] and OPAL [§] obtained data for cjl and ax separately and the latter 
collaboration even made a measurement of oa for the first time. The separation of ctl 
and ot is important because the former cross section enables us to extract the strong 
coupling constant a s and allows us to determine the gluon fragmentation density 
D g (x) with a much higher degree of accuracy as could be done before. Furthermore 
the measurement of a a provides us with information on hadronization effects since 
the QCD corrections are very small. 

As far as the theoretical achievements are concerned the order a s QCD corrections 
to the coefficient functions (k = T, L, A, £ = q,g), corresponding to the cross 
sections a k (0), have been calculated in the past in |10|, [□]] (see also 0). Also 
computed are the NLO corrections to the DGLAP timelike splitting functions Pij(x) 
= q,g) in | r2"f . Recently the order a 2 contributions to the coefficient functions 



Cl,i and Cr,i JT6| became available. From the latter one obtains the order 
a 2 corrections to the total longitudinal and transverse cross sections defined by 



s 



= (4) 



1 



which are equal to (k = T,L) 

Jo 



Cl q (z,Q 2 /M 2 ) + ±C k , 9 (z,Q 2 /M 2 ) 



(5) 



Here a^°'(Q 2 ) is the zeroth order annihilation cross section of process (0) which is 
identical to J2j<Joj in ec L- (2-14) of |J where / denotes a specific flavour (/ = 
u,d,s,c,b). Furthermore o" tot (e + e~ — ► X) = &t + &l which has been calculated up 
to order in | 17fl . In [|l^, |16[ it was shown that the order a 2 contributions to the 
coefficient functions are necessary to get agreement between the OPAL-data and the 
theoretical predictions. Since the OPAL-collaboration also measured the asymmetric 
cross section da^/dx in @) it will be of interest to compute the order a 2 corrections 
to this quantity too. In the QCD improved parton model it can be written as 

^M 2 ) = E jf t -W 2 ) ( D > (?0 - D " (?" 2 ) ) «"//•") 

(6) 

where / denotes the flavour of the quark. The asymmetry factor Af, which is defined 
in eq. (2.12) of ||, contains the products of the vector and axial vector electroweak 
couplings appearing in the 7 — Z and Z — Z interference term. It also includes 
the contribution of the photon and Z-boson propagators. The parton fragmentation 
densities denoted by D^(z,fi 2 ) depend in addition to the partonic scaling variable 
z also on the factorization scale //. Because of charge conjugation symmetry of the 
strong interactions we have used in (^) the identities 

fpNS _ ,rNS _ ,pNS ,pS _ ^NS _ n (7\ 

^-A,f — ~^A,f = ^A,qi A,q ~ ^A,g ~ U - \ 1 ) 

Up to order a 2 the non-singlet coefficient function C^ s g receives contributions from 
the following parton subprocesses where all quarks are taken to be massless 

O(a s ) : y-V+ff, (8) 
0(a s ) : V^"q"+q + g, (9) 
0(a 2 ) : V^"f +q + g + g, (10) 

V ->"<? +q + q + q. (11) 

Here the detected quark, which fragments into the hadron, is indicated by "g" and 
V — 7, Z. The above reactions also include the one- and two-loop corrections to 
process ([8]) and the one- loop corrections to process ([]). Further in reaction ([11|) the 
two anti-quarks can be identical as well as non-identical. In the case the anti-quark 
is detected one has to interchange q and q in eqs. (j^) - (JTTj) . 

The computation of the parton cross sections proceeds in the same way as has been 
done for the longitudinal Ci,i and transverse coefficient functions Cr,i presented in 
15] and |16| respectively. In the calculation one has to deal with the presence of ultra- 



violet (UV), infrared (IR) and collinear (C) divergences which have to be regularized 



2 



using the method of n-dimensional regularization. However there is one difference 
between the calculation of Cjy (k = T, L) on one hand and the computation of C^j 
on the other hand. This difference is due to the appearance of the 75-matrix in 
the interference term MyM\ + MaM v where My and Ma stand for the vector and 
axial-vector amplitude of the above processes. Here one has to find an n-dimensional 
extension for the 75-matrix occuring in Ma- For our calculation we have adopted the 
prescription for 75 given by 't Hooft and Veltman |H| (see also Breitenlohner and 
Maison ]20|). Since the axial vector vertex is represented by 7^75 one can simplify 



the traces using the identification 



6 



(12) 



which yields the same result as the prescription of 't Hooft and Veltman as is shown 
in |^T], P^j] . Although this prescription is consistent it has one drawback namely that 
the non-singlet axial vector current is renormalized in spite of the fact that it is 
conserved. Hence for each virtual correction where the 75-matrix appears in the loop 
one needs an additional renormalization constant to undo this unwanted effect. This 
constant has been calculated in ]24| and reads up to order a 2 



471"' 



C7 



5c 



+ 



47T 



Cl 



22 



CaCf { — 



44 107 



3- 



9 



(13) 



where the colour factors in QCD are given by Cf = (iV 2 — l)/2N, Ca = N, and 
Tf = 1/2 with N = 3 and the number of light flavours is denoted by rif. The 
rest of the calculation proceeds in the same way as performed for the deep inelastic 
coefficient functions C3 ^ j2B| which is the analogue of C^ s 9 . Apart from the check 



on the procedure outlined in [26] we will add a new one which has the advantage 
that we can get rid of the renormalization constant Za in ([13]) which is needed 
when the 75-matrix appears in the loop of the virtual Feynman graph. This check is 
based on the observation that the difference between the asymmetric and transverse 
parton cross sections does not contain the distributions denoted by 5(1 — z) and 
(ln fc (l - z)/(l - z))+. These singular functions originate from the one- and two- 
loop corrections to the Born-process (|]) and the contributions due to soft gluon and 
collinear fermion pair production in reactions (|9|)- (|TT|) . Hence these distributions 

Since the transverse coefficient function 

NS 



cancel in C^, Q 2 /u 2 ) - C^(z,Q 2 /fx 



C^ s q is known [JT6]] we can obtain C^ s g from the difference C^ s g — C^. The latter is 
only determined by the one-loop corrections to the regular part of process (hard 
gluon radiative part) and the regular part of (0), (0) (hard gluon radiation plus 
quark anti-quark production). Hence we only have to deal with the 75-matrix in the 
one-loop corrections to (§). However we have now checked that the following identity 
holds 





M^Mf + M^M^ + M^M^ 



3 



L v 



2M i A 1) *M^ ) 



+ 2M^ r M ( P 



(14) 



where Mjp is the order g l (a s = g 2 /An) contribution to the amplitude Mk (k = V,A). 

Here M^p and denote the amplitudes where the quark is attached to the vector- 
ial 

and axial- vector current respectively so that M A contains the 75-matrix. The term 
between the first pair of curly brackets on the left-hand side of (|bl| ) originates from 
the purely radiative process (§) whereas the term in the second pair of curly brackets 
refers to the interference between process (0) and the virtual corrections to (|9]). The 
latter is represented by the amplitude Mjp (k = V,A). Equation flU] ) reveals that 



one can get rid of the renormalization constant Za by shifting the 75-matrix from 
M A to the amplitude M A of the radiative process (|9|) so that this matrix becomes 
harmless. Actually one can now also choose the naive 75-prescription without altering 
the final result. The order a s corrections to C^r are calculated in 0] (see also eqs. 
(2.15), (2.16) in The order a 2 corrections calculated in this paper are presented 
as follows. First we split the coefficient functions Cj^ {k = T, A) in two parts i.e. 



T,q > ^T,q > 



(15) 



,pNS 
^A,q 



NS,id 



CNS,nid 
A,q ~ ^A,q 



(16) 



where the last term in the above equations is only due to the contribution from 
identical anti-quarks in reaction The second order contributions to both parts 

see appendix in A in [TE\) as follows 



can be now obtained from C^' md and C^' ld 



,pNS,nid,(2) 



f pNS,nid,(2) 
^T,q 



4vJ 



2(1 - z){2 In z - 4 ln(l - z) + 1) 



O 2 

In + 4(1 - z)(AS 1 2 (1 - z) - 8Li 3 (-z) - 4((2) ln(l - z) + 4 lnzLi 2 (-z) 
fi 2 



25 



24 8 



+ 3Li 2 (l - z) - W{1 -z)- In zln(l - z) + — ln(l - z)) + (— + 

2 5^ z 



24 



16 - 16z + 8z 2 + — z d )(Li 2 (-z) + lnzln(l + z)) + (-24 + 8z + 8z 2 



24 , N> , s , , 12.,, . 24 2 202 24 

+ -z 3 ( 2 + 10 -2z- Az 2 - -z 3 ) In 2 z + -— ^ + - + — z - - 

5 5 bz z 5 5 5 



24 9 9 24 2 



22 O 2 

+ C A C F { — (1 - 2) In + 4(1 - z)(4Li 3 (-s) - 2£ 1)2 (1 - z) 



A 



25 12 
+ 2C(2) ln(l - 2) - 2 In zLi 2 (-z) - — ln(l - z)) + (- 



6 



5z 2 



- + 8 + 8z 



4z 



2 ^ z 3 )(Li 2 (-z) + In z ln(l + z)) + (4 + 4z - 4z 2 - ^z 3 )C(2) 
o 



■? 6 2 ,12 182 398 12 2 . , 12 

+ -2 - 2z + 2z 2 + -z 3 ) In 2 z + — + — —2 + — z 2 ) lnz - — 

5 5z 15 15 5 5z 



823 823 



12 



45 + iT + T" 



8 O 2 8 19 

+ n/OT, <{ - -(1 - z) ln^ - -(1 - z)(ln(l - z) + lnz - 



6 



(17) 



(p NS,id,(2) _ ( pNS,id,(2) 
A,q ^T,q 



8{l + z){AS li2 {-z)-2U 3 



+ 4 ln(l + z)Li 2 (-z) + 2C(2) ln(l + z) + 2 In z ln 2 (l + z) - In 2 z ln(l + z) 



-2<(3)) + (£-$-« 



24 

2 1 — z 3 )(Li 2 (-z) + lnzln(l + z)) + (8-8z 
o 



24 1 2 24 7? 

- 8z 2 + #z 3 )C(2) + (-4 + 4z + 4z 2 - ^z 3 ) In 2 z + (-^ + ^ 
5 5 5z 5 



72 



24 



+ — z -z 2 ) lnz + — + 

5 5 5z 5 



24 104 104 



24 



5 5 

where the Riemann zeta-function ((n) and the polylogarithms Li n (z), S njP (z) can be 
found in P7J. 

Notice that in C^' ld (k = T, A) we have omitted contributions as represented by the 
cut graphs in fig. ]l|. The photon cannot couple to the cut fermion triangle because 
of charge conjugation invariance. However the Z-boson decouples too if one sums 
over all flavours in one family. This is because the Z is connected to the quarks 
via the axial-vector coupling constant representing the weak isospin component 
of a specific flavour / with the property X^/=u = ®- From now on we will 
assume that in the inclusive state one sums over all members in one family so that 
the contribution due to fig. [I] can be dropped. 

The first quantity we would like to study is the flavour asymmetry sum rule which is 
defined in eq. (2.23) of 0. It is given by 



Z Q A = J2MQ 2 ) [ 1 dz 1 Q ( P ( Df{z u f)-Df{z ut ? 

H,f J ° V 



■ t dz 2 C^ q {z 2 ,Q 2 /^) 
Jo 



(19) 



( f) 

where is a conserved additive quantity. The first moment of the non-singlet 
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coefficient function calculated up to order a 2 is equal to 

•1 2 



/ dz 2 Cf q (z 2 ,Q 2 / f , 2 ) = l 
Jo 



a s (Q z 
An 



12/5 C F C(3) 



+ 



a s (Q 2 

47T 



where (3q is the lowest order coefficient of the beta-function given by 



~ 3 T f n f . 



,(3) 
-A,q 



(20) 



(21) 



Notice that the first moments of q and Dj? — D¥ are separately scheme indepen- 
dent. Further there are no order a s corrections to the first moment of C^ s g || and 
the order a 2 correction is proportional to (3q. A comparison between (|20|) and the 



order az corrected R, 



ee — °toH e e 



X) I <j(°> [[IT] reveals that the coefficients of the 
Riemann zeta-f unctions (here ((3) only) are exactly the same. Furthermore if one 
drops all rational numbers in R ee one obtains exactly (f20|). Following the arguments 
in |30| one can make an interesting conjecture about the third order term c A which 



has not been calculated yet. Suppose that all rational numbers in c A ' q are zero and 



that the coefficients of the Riemann zeta-functions C( n ) (here £(3) and C(5)) are the 
same as in R ee then we can make the following conjecture 



.(3) 
"A,q 



C A C 



+ C 2 F T f n f 



+ C F T]n) 



572C(3) + 880C(5) 



304C(3) - 320C(5) 



C F C\ 



10948 440 
^p(3) - ^C(5) 



+ CAC F TfUf 



7168 w , 160 w , 

-g-C(3) + — C(5) 



1216 
~9~ 



C(3) 



N 



K(3) 



'.22) 



where d a b c denote the structure constants which emerge from the anti commutation 
relations of the generators of the group SU(N). We now want to study the effect of 
the order a 2 correction on the asymmetric fragmentation function and compare the 
result with the OPAL data ||. The fragmentation functions F% will be defined by 
(see FT) 



F*(x,Cp) 



1 d<jf;(x, Q 



(k — T,L, A). 



o~tot dx 

If we sum over all hadrons of species H we obtain the quantities 

F k (x, Q 2 ) = Fkix, Q 2 ), (k — L,T), 

H 

F a (x,Q 2 ) = J2QhF^(x,Q 2 ), 



(23) 

(24) 
(25) 



H 
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where the sum in ( pq) is taken over all charged hadrons. From (ra) and ([24]) we infer 
that Fa gets only contributions from the valence fragmentation densities Dy } = 
— . Hence the measurement of Fa provides us with information about the 
x-behaviour of the valence fragmentation functions. The latter are available for H = 
7r , K ,p,p in JHJ where they are parametrized in leading log (LL) and in next-to 



leading log (NLL, MS-scheme). For H 



71 



K + , p we obtain from (|6|) and 



F%(x,Q' 



Otot 



1 dz 

z 



V Z / V Z 



(26) 

with U = u (tt + ,K + ) and D = d (7r + ) or D = s {K + ). The proton contribution is 
given by F% = 0.16 FJ + where the factor 0.16 originates from [32] where one has 
estimated Fp+p = F[ +p + F T +P by putting Fp+p = 0.16 F 7r++7r " . Since in gg one 

26f) is negative over the whole x- 



has taken Dy V 



Dy D we observe that Ff (x, M%) 



region. This property can be traced back to the value of the electroweak angle leading 
to Ad(M§) / Ajj(Mz) ~ 2. Therefore all hadrons with positive charge {Qh > 0) give 
a negative contribution to Fa{x, M§) (p5|) . If is the anti-particle of H we have the 



relation F^ 



-F H 
r a ■ 



Because of Qg = —Qh in Q25|) the anti-particles , K~ ,p) 
also give a negative contribution to Fa(x, Q 2 ). Therefore the parametrization in p2| 
predicts a negative Fa(x, Q 2 ) (|25| ) over the whole x-region at Q 2 = M|. 
In our plots discussed below a comparison will be made with the OPAL data || so 



Q 2 in (|5D and n/ = 5. 



that we have to choose Q 2 = M|. Further we take /i 2 

The running coupling constant is chosen to be a; a (M|) = 0.126. Finally we want to 
emphasize that a full next-to-next-to-leading (NNLO) analysis of Ft and i 7 ^ is not 
possible yet because of the missing three-loop contributions to the DGLAP splitting 
functions. Therefore the order a 2 correction, which can be only attributed to the 
coefficient functions in 



'S 

and 



16 



have to be considered as an estimate. The 
order a 2 corrected Fl is complete because here the NLL fragmentation densities and 
the order a 2 



m). 



corrected coefficient functions are available (see 
In fig. | we have plotted F%°, F^ LO and Ff NLO together with the OPAL data 
(see also fig. 4 in ||). There is a difference between F^° and F^ LO but the order 
a 2 corrections shown by F^nlo are uno bservable. Furthermore the theoretical curves 



are above the data, 
for the ratio 

Ra(x,Q 2 



In fig. 8 of || the OPAL-collaboration also presented the data 



F A (x,Q> 
F(x,Q 2 



F(x,Q 2 ) =F T (x,Q 2 ) = F L (x,Q 2 ). 



(27) 



In fig. | these data are compared with the theoretical predictions Ra°, Ra L °i an d 
R^ NL °. Here we see the same features as has been observed for Fa in fig. 0. There 

NNLO 



is no difference between Ra LO and Ra NLU and only the order a s corrections, rep- 
resented by Ra LO , are visible. Also in this case the data are below the theoretical 
predictions. 

From the data one can infer the integrated fragmentation function for which the 
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theoretical predictions corrected up to order a 2 are given below 



i 



dxF% NLO {x, M\) = -0.016 (-0.023), (28) 



o.i 



,i i 

/ dx-xF% LO (x,Ml) = -0.0020 (-0.0027). (29) 
io.i 2 

The experimental values for fl28|) and © are -0.0229± 0.0044 and -0.00369± 0.00046 



respectively. Since the fragmentation densities in have a limited range of validity 
we have imposed a lower bound on the integration which is given by a; = 0.1. Between 
the brackets in fl28|) , fl29f ) we have quoted the LO results. It turns out that the 
latter are in better agreement with experiment than the NLO and NNLO numbers. 
Further the values of the integrals also hold in NLO since the order a 2 corrections 
are extremely small. 

The OPAL-data indicate that at low x Fa{x,M^) might become positive. If this is 
the case one has to assume that in this region Dy V (x, jj 2 ) > Dy D (x, fi 2 ) provided the 
zeroth order contribution to C^ s ? which is given by 5(1 — z) dominates the integral. 
Summarizing the above we conclude that the order a 2 corrections to Fa are negligible 
and we do not expect that this will change when the effect of the three-loop DGLAP 
splitting functions are taken into account. Furthermore the above results reveal that 
the measurement of Fa puts some constraints on the valence fragmentation densities. 



In particular it means that the NLL parametrizations in |32| have to be modified in 
order to get agreement with the OPAL data. 
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1 Figure captions 



Fig. 1 Diagrams with four quarks in the final state containing a cut triangular quark 
loop. 

Fig. 2 Contributions to the asymmetry fragmentation function F^{x, Q 2 ) ( p5f) at 



Q = Mz using the fragmentation density set of [32|. Solid line: LO. Dashed line: 



NLO. Dotted line: NNLO. The exprimental data are taken from OPAL ||. 

Fig. 3 The ratio Ra(x, Q 2 ) p7| ) at Q = M z using the fragmentation density set of 
|32j. Short dashed line: LO. Solid line: NLO. Dotted line: NNLO. The exprimental 



data are taken from OPAL 
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